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ABSTRACT 


The  effects  of  a  neutron  charge  distribution  are  considered. 
It  is  suggested  that  the  charge  distribution  is  of  the  form 

The  effects  of  this  charge  distribution  for  various  values 
of  the  parameters  c  and  f,  on  the  mean  square  radius  of  carbon, 
and  on  the  charge  form  factor  of  carbon  at  high  q  are  shown.  The 

effect  of  the  neutron  becomes  appreciable  in  carbon  above  the 

2 
diffraction  minimum.  The  neutron  in  this  range  of  q  adds  approx- 
imately 50%  to  the  carbon  form  factor. 
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1.   INTRODUCTION 

A  review  of  current  literature  concerning  the  nuclear  charge 
distribution  shows  that  the  effect  of  the  neutron  is  largely  ignored 
except  in  deuteron  studies.  This  paper  proposes  to  suggest  what 
effect,  if  any,  a  finite  neutron  will  have  on  the  distribution  of 
charge  in  any  nucleus. 

Before  the  charge  distribution  in  a  nucleus  can  be  calculated, 
the  wave  functions  describing  the  motion  of  the  nucleons  and  the 
nucleons'  internal  charge  distributions  must  be  known.  The  proton 
charge  and  magnetic  moment  distributions  and  the  neutron  magnetic 
moment  distribution  are  known.  The  neutron  charge  distribution  is 
not  well  known,  however. 

If  wave  functions  are  known,  then  the  neutron  charge  distribution 
may  be  related  to  the  measured  nuclear  charge  distribution.  In  this 

paper  the  diffraction  minimum  observed  in  electron  scattering  by  a 

12 
C   is  used  to  determine  the  nucleon  wave  functions. 


2.   THE  NUCLEAR  FORM  FACTOR 


In  order  to  develop  a  theory  for  the  charge  distribution 
of  a  nucleus,  it  is  assumed  that  each  nucleon  maintains  its 
identity  in  the  nucleus.  This  assumption  leads  to  models  which 
successfully  predict  various  experimental  results.  The  assumption 
implies  that  each  nucleon  in  the  nucleus  may  be  described  independ- 
ently of  the  other  nucleons.  The  nuclear  wave  functions  may  be 
written  as  follows: 

f'lt)  V\l) fit 


(2 


-«  f(t--/0=  A 


Ai 


iKCr.) 


Hiti) 


t'tf) 


r.  describes  the  position  of  the  center  of  mass  of  the 


.th 


i   particle.  The  ^'  implicitly  include  a  spatial  wave  function, 
hereafter  referred  to  as  ^ /  spin  and  isospin  functions .  All  of 
these  functions  are  orthonormal.  The  use  of  the  othonormality 
will  be  implicit  in  what  follows.  The  total  wave  function,  x  , 
is  also  normalized  and  is  antisymmetric  on  interchange  of  any  two 
nucleons  (a  consequence  of  the  nucleons  being  fermions) . 

The  purpose  of  this  paper  is  to  examine  the  nuclear  charge 
distribution  a (R) .  This  cannot  be  measured  directly.  What  can 
be  measured  is  the  differential  cross  section  for  scattering  electrons, 

dcr 

dSL' 


The  cross  section  may  be  separated  into  a  term  due  to  the 
charge  distribution,  and  a  term  due  to  the  magnetic  moment  of  the 

nucleus  [1] .  The  charge  distribution  term  leads  to  the  charge  form 

2  ... 

factor,  F(q  ),  which,  in  Born  approximation,  is  the  Fourier  trans- 
form of  the  normalized  charge  distribution. 


(2-2)  F(%%)  =  -^  fe^'W?)c{3r 


The  Born  approximation  will  be  used  throughout  this  paper. 


The  nuclear  charge  distribution ,/0  (R) ,  is  determined  by  the 
nuclear  wave  functions  and  the  charge  distribution  of  the  nucleus. 


1* 

r 


(2-3)  f(*>  -.  f^e^-x*)/'  jprt,-  A)d 

where  0  is  the  nucleon  charge  operator: 

(2-4)   2  r  £  /  <  &  -  Vy , 
/       i.-i  ^ 

The  f.  are  the  charge  operators  for  the  individual  nucleons. 
Putting  equations  (2-3)  and  (2-4)  into  equation  (2-2)  yields 


(2-5)  rv)«  -jr 


e 


Let  <k  (r . )  be  denoted  by  j.  . .  Note  that  after  integrating  over  all 


'ID 

■hVi 

except  the  i       coordinate,  we  have 


(2-e,    ma-^ifr    :    i.^/j^.A    . 


We  see  from  equation  (2-5)  and  (2-6)  that,   if  the  single 
particle  wave  functions  are  orthogonal,  the  charge  distributions 
for  all  of  the  protons  are  equal,  those  of  all  the  neutrons  are 
equal  and  the  charge  distributions  are  additive. 

Using  (2-6)  in  (2-5)  implies 

The  relation  between  R  and  r.  is  R  =  r .  +  r .  • 

i        11 

Insert  the  following  equation  into  (2-7) : 


(2-8,//(C- (*-*))  A' 


Equation  (2-7)  becomes 


1  tl\    .i\ 


(2-9) 


F(f) ..  4-  fatj^WiW-W& *  <• 


(2-10) 


Integrate  over  R, 


i  . 


Separate  integrals  into  one  over  r.  and  one  over  r. 


.U](Aii%^v)^,(r^ni^)idVl)j 
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Split  the  sum  into  one  over  protons,  f .  =  f  .and  one  over  neutrons, 

1    p  ' 


f.  =  f 
1    n 


FO-- 


Define^F  =    (      l\   V  f    ^  >N  ,  } 


(2-12) 


U*  /*tVV»c?'M 


r 


,2-n,  ijbKv^<?HjV. 


Since  we  have  integrated  over  r,  F. .  =  F.  for  all  i.  Use  (2-12) 
and  (2-13)  in  equation  (2-11) . 


(2-14) 


.   F  -'-i:  F  +  A.  F  J-  £  F 


Define 


(2-15)   P 


Cr*-*~         A         j., 


Then  write  equation  (2-14)  as 
(2-16) 


Rf)  ■  [  >>  +  i  S]F* 


*«i4 


The  nuclear  form  factor  has  been  separated  into  three  parts: 
(1)  F  ,  the  form  factor  of  the  proton 


11 


(2)  F  ,  the  form  factor  of  the  neutron,  and 

n 

(3)  F   ,.,  a  form  factor  due  to  the  distribution  of  the  centers  of 

cm-A 

mass  of  the  nucleons. 

F  is  valid  for  any  nucleon  distribution,^  (R) ;  however,  >  (R) 
may  be  expanded  (since  any  function  of  angle  may  be  expanded  in  terms 
of  spherical  harmonics)  as  follows: 


<2-i8,  ^e^  --    f^^y^) 


(  ft  «ftt,/*=o 

where  K  i.sf   (R)  averaged  over  all  angles  or  the  spherically 

symmetric  part  of  /^(R) . 

2 

F(q  )  may  be  expanded  in  terms  of  the  multipole  expansion  of  the 

charge  distribution.  This  yields  an  electric  monopole  term  plus  higher 
order  multipole  terms  in  the  form  factor.  The  simplest  nuclei  to 
analyze  are  then  the  spherically  symmetric  nuclei,  since  these  have 
only  the  electric  monopole  term.  All  further  work  in  this  paper  will 
deal  only  with  this  type  of  nucleus. 


12 


3.   MEAN  SQUARE  RADIUS  DETERMINATION 

A  quantity  which  is  experimentally  easy  to  determine  is  the 
mean  square  radius.  We  shew  below  that,  theoretically,  this  radius 
may  be  expressed  in  terms  of  the  particle  wave  functions  and  in  terms 
of  the  charge  distributions  of  these  particles. 

The  form  factor  for  a  spherically  symmetric  charge  distribution 
4(r)   is  ^  ^ 

(3-D    R%T)=  fcl*s"ry>  ^r 

©  f 

We  have  on  expanding  sin  (qr) : 
This  implies  that 

J  Ftp  \   ,   .-  -  <-£?    :S> 

(3-3)    J  ^T     'tra         ^ 
From  equation  (2-16)  we  have 

,2-i6,    Fee1)  -  [  f,  +  £  5  ]  £-   • 
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Using  this  in  equation  (3-3)  yields 

(3-4)  <r 


^  su 


r 


5  O 


In  equation  (3-4)  we  have  used 

The  directly  measurable  quantities  in  equation  (3-4)  are 
the  mean  square  radii  of  the  nucleus  and  of  the  proton.  The  mean 
square  center  of  mass  radius  is  dependent  on  A  and  the  configuration 
of  the  nucleons  in  the  nucleus.  This  configuration  may  be  different 
for  different  elements  of  equal  A. 

It  is  thus  possible  to  write  the  nuclear  form  factor  in  terms 
of  the  nucleon  form  factors  and  nucleon  wave  functions  which  describe 
the  states  of  the  nucleons.  See  equation  (2-16).  However,  starting 
with  experimental  data  on  the  nuclear  form  factor,  it  is  impossible 
to  separate  the  contributions  due  to  the  nucleon  charge  distributions 
and  that  due  to  the  nucleon  motion  unless  a  model  for  the  nuclear  or 
nucleon  structure  is  assumed.  That  is,  the  analysis  will  be  model 
dependent.  The  validity  of  any  model  is  measured  by  how  well  it 
agrees  with  experiment. 
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4.   POSSIBLE  NUCLEON  CHARGE  DISTRIBUTIONS 


There  has  been  a  great  deal  of  work  done  on  the  proton 

which  suggests  that  it  has  a  gaussian  charge  distribution  [2] . 

.1 


(4-1)   I  ( r)    *   "-J -   C 


-* 


The  form  factor  corresponding  to  this  charge  distribution  is 

(4-2)  rep  -  t. 


Figure  1  shows  a  comparison  of  this  form  factor  to  experimental 
data  [2] .  We  assume  here  that  f  is  correct,  and  that  the  charge 
distribution  and  the  matter  distribution  are  identical.  With  this 
assumption  it  seems  valid  that  the  matter  distribution  in  the  neutron 
could  also  be  distributed  in  the  same  manner.  That  is  to  say,  the 
matter  distribution  of  the  neutron  is  a  gaussian  with  approximately 
the  same  radius  as  the  proton.  It  is  knewn  that  the  total  charge 
of  the  neutron  is  zero.  However,  its  magnetic  moment  suggests  that 
the  neutron  does  have  a  charge  distribution.  By  analogy  with  the 
proton  the  following  form  is  appealing: 


(4-3) 


(  i  /o  c%  C 


This  is  a  double  gaussian,  each  with  unit  charge.  The  negative  charge 
has  the  larger  radius. 

Our  assumption  that  the  matter  distribution  of  the  proton  is  the 
same  as  that  of  the  neutron  implies  that  the  radius  of  the  neutron's 
negative  (larger)  gaussian  is  approximately  that  of  the  proton. 
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f+A 


fr-e 


o<L 


s  Ji^ 


X    -j>**e1<    J3kM,     ? /;.     -'si«s 


*   . 


X    \ 


_     I 


I  • 


0        ' 
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Figure  1 
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We  take  f  =  .44F  (slightly  larger  than  the  proton) .  This  leaves 
one  parameter  to  fix  in  the  neutron.  The  free  parameter  will  be 
fixed  in  this  paper  by  using  the  harmonoc  oscillator  model  and 
experimental  carbon  data.   (For  alternative  method  see  Appendix  C.) 
The  form  factor  correspond  to  f  is,  in  Born  approximation 


(4-4) 


# 


Figure  2  is  F  for  various  values  of  c  compared  to  experimental 

data  [3].  This  data  does  not  deny  any  of  the  suggested  curves,  but 

favors  F  =0. 
n 
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e     <S 


Figure  2 
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5.   HARMONIC  OSCILLATOR  WAVE  FUNCTIONS 

2 

In  order  to  determine  <  r  P    we  are  forced  to  assume  a  model 

n 

to  predictor  f   .  The  harmonic  oscillator  potential  is  a  simple 
potential  which  appears  to  agree  well  with  experiment.  This  theory 
has  the  added  advantage  that  its  parameters  are  fixed  independently 
of  the  mean  square  radius  of  the  charge  distribution. 

Consider  non-interacting  particles  in  a  spherically  symmetric 
harmonic  oscillator  potential  well: 


\/tr)    --  -7-  ^^X^ 


(5-1)   V  I  r  j  ~-  y 

The  stationary  states  or  orbitals  occupied  by  the  particles  are 
solutions  of  Schrodinger ' s  equation. 


(5-2) 


4   yn 


Using  equation  (5-1) ,  and  defining  a  =  *  a  ,  equation  (5-2)  can  be 

m  w 
written  as 


JL  iJ  [~  vl  + 

(5-3)     7         »v\      L  ^ 


£}]^(r)  =  Ef(?) 


(5-4)   where  v  \  * 


?l 


Assume 


(5-5) 

Using  equations  (5-4)  and  (5-5)  in  equation  (5-3) 
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(5-6)  R  or  ^l  &r   J-1   -v  \       ' 

l     [     '        jL.  f  3,^iZ  \  ,     L    iLi  1   = 

The  solutions  to  the  angular  equation  are  the  spherical  harmonics. 
(See  Appendix  B.) 

For  convenience  in  solving  the  radial  equation,  define 

(5-7)   "R(v)  s   ^Jj[} 

V 

The  radial  portion  of  equation  (5-6)  may  be  written  as 

(5-8)  eTucO 


The  solution  of  equation  (5-8)  and  equation  (5-7)  yields:   (See 
Appendix  B.)  ,1 


,s-9)    R^.r'e."  ZC<(*)  *tK. 

These  equations  also  yield: 
(5-10)   ^  =  f^J  [^|f  f  til 

Where  A  =  2N  +  1  denotes  the  energy  levels.  The  harmonic  wave 
functions  are  then  (See  Appendix  B.) 

(5-n)  t<,r  R*X.  . 
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The  above  solutions  are  based  on  the  central  potential  having 
a  fixed  origin.  In  reality,  the  central  potential  should  be  measured 
from  the  mean  position  of  the  nucleons.  It  has  been  shown  that  this 
complication  introduces  a  factor  £  *'*   into  the  form  factor  [4]  . 

From  equation  (2-4)  and  the  above  consideration  we  have  (See 
Appendix  B . ) 

,5-xa,     f^=   e^/V'|V. \cir 

From  equation  (2-16)  we  have 
(5-13)   F~     s  —  P~    p 

If  the  harmonic  oscillator  model  is  applicable  and  if  the 

nucleon  form  factors  F  and  F  and  the  nuclear  configuration  &•©•# 

n     p  r 

n,  1,  m  must  be  known  for  each  nucleon)  are  known,  then  the  total 
nuclear  form  factor  may  be  written.  For  example,  for  p-shell 

nuclei  (4<:A£16,  2^  Z^  8)  [5]: 

A 


(5-14) 


j_  (a  -  (a-<)  4T  )c 

-  A   v  y 


The  F  .   from  Appendix  B  have  been  used.  In  Carbon  A=12,  Z  =  N  =  ^- 
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Or  the  total  form  factor  is  ,      g.1  -t  /  i  .  Ar    \ 

shv(f,*eX«- V;* 

0  rr 

F    has  a  zero  at  q    =-*»•.     "Therefore,  our  cross  section  measurements 
c  z 

a 

should  also  show  a  zero  at  this  point.  Experimental  data  show  a 

2 

marked  minimum  in  the  cross  section  which  can  be  used  to  fix  a  . 
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6.   COMPARISON  TO  EXPERIMENTAL  DATA 

In  the  harmonic  oscillator  model,  the  center  of  mass  form 

factor,  F    ,  is  completely  determined  by  the  position  of  the 

2  2      2 

minimum  in  a  plot  of  F  (q  )  vs.  q  . 

Carbon  offers  a  very  good  method  of  fixing  the  neutron  parameters 
since  it  is  spherically  symmetric  and  the  experimental  data  show  a 


prominent  diffraction  minimum  [6] . 


,  ,x    -&\.-V) 


(6-1) 


(6-2)   A   -   "TT-     m  '   3.Z"?  F 


Figure  3  is  a  comparison  of  equation  (6-1)  to  experimental  values 
for  the  neutron  parameters  of  Figure  2  [6] .  This  figure  indicates 
that  the  best  neutron  form  factor  is  given  by 


- 1^'      -  i'^D 

(6-3)    K  -   C         —  C 


(6-4) 


The  mean  square  radius  of  carbon  is  then 

a  .6"  F' -  ,U  Fl  +  5,41  F'«  Cot  F' 

This  radius  is,  of  course,  model- dependent  because  of  the  manner 

in  which <r  ^_  has  been  found.  However,  it  agrees  well  with  current 
cm-A  r 

model  independent  determinations  of  ^  r  'c  at  Darmstadt  in  which 
^r2^c  =  .587  +  .19  F2  [7]. 


23 


CD  £=o    v 
CpF^.saT) 


Figure  3 
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7 .   DISCUSSION 

The  preceeding  results  suggest  two  alternative  approaches 
for  determining  the  validity  of  the  neutron  radius  in  equation 
(6-4) .  First,  assume  that  the  nucleons  are  described  by  harmonic 

oscillator  wave  functions,  in  which  the  range  parameter  a  is  adjusted 

2 

to  give  the  experimentally  observed  diffraction  minimum  in  F  (q  )  . 

This  allows  us  to  write 

2 

We  would  then  use  the  model  independent  <rx?^  determined  at  low  q  , 

the  accepted  value  of  ^  r  ?  determined  from  IF"  data ,  and  the  shell 

model  prediction  of  <r  r      .  Extreme  accuracy  in  the  measured 

2 

quantities  is  required  due  to  the  small  size  of <  r  >  . 

The  second  alternative  would  be  to  assume  that  the  nucleon 
center  of  mass  distribution  is  the  same  for  elements  with  equal  A 
but  different  Z.  This  assumption  is  based  on  the  charge  independence 
of  nuclear  forces.  The  assumption  could  be  checked  in  a  model-' 
dependent  way  by  seeing  whether  or  not  their  diffraction  minima 
occurred  at  the  same  value  of  q.  Measuring  the  model- independent 
radii  of  two  elements  with  the  same  atomic  weight  at  low  q  and 
subtracting  the  following:  ~ 


(7-3) 
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Thus  C  r  S  is  obtained  in  terms  of  only  the  model  independent 
mean  square  radii.  The  assumption,  however,  that^  r  ^a  =<^r  _£ 
would  only  be  valid  if  the  nucleon  configurations  of  the  elements 
were  identical.  One  possibility  would  be  to  use  chromium-50  and 
titanium-50.  These  two  nuclei  have  total  angular  momentum  (J)  equal 
to  zero.  Therefore,  they  have  only  a  monopole  form  factor.  Any 
other  pair  of  mirror  nuclei  would  do  equally  well. 

Both  of  these  methods  could  be  checked  by  comparing  the 

theoretical  form  factors  to  the  experimental  form  factors  at  high 

2 

q  .  This  is  the  manner  that  the  neutron  parameters  were  fixed  in 

this  paper. 

The  limitation  on  the  parameter  f  in  the  neutron  charge 
distribution  (that  the  neutron  matter  distribution  equals  the 
negative  charge  distribution)  is  of  questionable  validity.  If 
this  restriction  is  removed,  and  we  use  instead  the  electron- 
neutron  interaction  condition, 

\r  \ 

(7-5)  %JJL  \    ,   r  ,02.1   r       [8],  [9], 
which  implies 

(7-6)     <  r*  >    =  -  JH  Fl 


2 

Then  the  difference  between  the  two  neutron  parameters  is  .084  F  . 
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F„=£ 


When  this  restriction  is  used  and  the  form  factor 


is  again  fitted  to  various 


2       2 

experimental  data,  a  good  fit  is  obtained  for  c  =  .09  F  , 

2        2 
f  =  .174  Fz  . 

This  choice  of  parameters  leads  to  a  form  factor  which  gives 

2         2  -2 

a  smaller  radius  at  low  q  ,  but  at  q  of  about  25  F   the  radius 

is  much  greater  than  with  the  previous  choice  of  parameters.  The 

form  factor  of  the  neutron,  with  this  choice  of  parameters,  becomes 

2      -2 

about  three  times  that  of  the  proton  at  q  =  25  F  .  This  doesn't 

seem  reasonable.  The  previous  selection  of  neutron  parameters  is, 
consequently,  preferable. 
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8 .   CONCLUSION 

Current  experimental  data  point  to  the  possibility  that  the 
charge  distribution  of  the  neutron  has  a  greater  mean  square  radius 
than  previously  thought.  It  appears  that  a  value  of<rn^=  -.2F 
is  not  unlikely.  In  fact,  such  a  value  is  supported  by  current 
measurements  of  nuclear  radii. 

It  should  be  noted  that  it  is  presently  impossible  to  separate 
the  effect  of  the  neutron  from  the  effect  of  the  wave  functions  in 
currently  available  experimental  data  on  electron  scattering.  A 
modification  of  the  wave  functions  can  give  form  factors  which  fit 
the  data  over  a  large  range  of  q  without  considering  any  neutron 
effect. 

The  approach  suggested  in  this  paper  has  the  advantage  of 

simplicity.  Simple  wave  functions  and  a  siitple  neutron  give  agree- 

2 

ment  with  F(q  )  data. 

If  the  charge  distribution  of  the  neutron  can  be  determined 
independently  (for  example,  through  neutron-electron  interactions) , 
then  this  information  could  be  used  in  fitting  wave  functions  to 
experimental  data. 


28 


APPENDIX  A.   COMMON  INTEGRALS 
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APPENDIX  B 


HARMONIC  OSCILLATOR  WAVE  FUNCTIONS  AND  FORM  FACTORS 


1.  Spherical  Harmonics 
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2.  Radial  Wave  Functions 
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3.  Harmonic  Oscillator  Wave  Functions    t  ^* 
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4.  Electric  Monopole  Form  Factors  for  Independent  Particles 
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APPENDIX  C 
CLASSICAL  MAGNETIC  MOMENT  DETERMINATION 
Classically  the  magnetic  moment  m  is  defined  as 

where  r  is  the  distance  from  some  point  of  interest  and  J  is  the 
current  density.  For  an  arbitrary  charge  distribution,  /^(r),  rotating 

A 

about  an  axis  with  velocity  CJ  K 


(c-2)   J  ^  C  ^  *  r  )  f?  ir 


Using  this  in  equation  (C-l)  yields 


<c"3>  *  *iz; f* rVs  * ?r*c?) d  ir. 

Using  the  vector  identity   "f  K  ^  K  y  -  (y    T  )  UJ  ..  (  r  w  )  y^ 
this  becomes 

(c-4)  ^  =  J-  f(   r*  u  _  y»*u  £  )>£?)  r|V 

We  are  interested  in  finding  the  magnetic  moment  for  neutrons 
and  protons  projected  onto  the  R  axis.  For  a  spherically  symmetric 
charge  distribution  we  have 

(C-5)   yyv,  K  -s   

?  C 


\Yl<+>j*t*)c\  >'  -  j  r?<->  to±*si>(f)c\*r 


The  second  term  integrates  to  zero  when  the  angular  integration  is 
performed. 
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So  the  magnetic  moment  about  the  Z  axis  is  for  the  proton 

(c-7)   M   =  J^t    <l.  r^  "? 

and  for  the  neutron 

(0-8)        |fi,%  =       t?L     £   Y?  > 

The  ratio  of  the  proton  moment  to  the  neutron  moment  is  then 

(C-9)    IS   ~    — 


assuming   £j  -  U)a 

2  -» 

This  agrees  in  sign  with  the  value  predicted  for^  r  f   as 

2 

compared  to<"rp>but  is  only  the  correct  order  of  magnitude  and, 

therefore,  gives  neutron  parameters  larger  than  those  in  this 
paper. 

(do,  <v;>  =  --£.    =  -----   -  -.«^*. 
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